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Abstract 

We argue that quantum-gravitational fluctuations in the space-time background give 
the vacuum non-trivial optical properties that include diffusion and consequent uncertain- 
ties in the arrival times of photons, causing stochastic fluctuations in the velocity of light in 
vacuo. Our proposal is motivated within a Liouville string formulation of quantum grav- 
ity that also suggests a frequency-dependent refractive index of the particle vacuum. We 
construct an explicit realization by treating photon propagation through quantum excita- 
tions of L'-brane fluctuations in the space-time foam. These are described by higher-genus 
string effects, that lead to stochastic fluctuations in couplings, and hence in the velocity of 
light. We discuss the possibilities of constraining or measuring photon diffusion in vacuo 
via 7-ray observations of distant astrophysical sources. 
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1 Introduction 



The propagation of light through media with non-trivial optical properties is subject 
to three important effects: one is a variation in the light velocity with photon energy, 
namely a frequency- dependent refractive index, a second is a difference between the 
velocities of light with different polarizations, namely birefringence, and the third is 
a diffusive spread in the apparent velocity of light. We have argued previously that 
quantum-gravitational fluctuations in the space-time background may endow the 
conventional particle vacuum with such non-trivial optical properties, in particular 
a frequency-dependent refractive index We have also observed that such an 
effect may be severely constrained by careful observations of distant astrophysical 
objects whose emissions exhibit short time structures, such as Gamma-Ray Bursters 
(GRBs) 1^. The possibility of quantum-gravitational birefringence has been raised 
within a loop approach to quantum gravity 0. The purpose of this paper is to 
propose how quantum-gravitational diffusion may spread the arrival times of photons 
from distant sources, even if they have the same energies (frequencies), corresponding 
to stochastic fluctuations in the apparent velocity of light. 

An example of this phenomenon in a conventional optical medium, which illus- 
trates clearly the intuition behind our proposal, has been discussed in Light 
propagating through ice may encounter air bubbles, which have a different refrac- 
tive index and hence induce scattering and diffusion. We have argued previously |l|] 
that foamy fluctuations in space-time generate a quantum-gravitational 'aerogel' 
with an effective refractive index ~ E/Mqg in vacuo, where Mqg is some char- 
acteristic scale that may approach the Planck mass Mp. Now we argue that the 
quantum uncertainties in these foamy quantum-gravitational fluctuations in turn in- 
duce fluctuations in the refractive index that act as scattering centres analogous to 
the air bubbles of 0. As illustrations of these quantum-gravitational fluctuations, 
one might consider foaming microscopic black holes that induce local fluctuations in 
the wave front close to their horizons, leading on larger distance scales to diffusive 
broadening of any light pulse. 

In this paper, we provide a mathematical formulation of this physical intuition, 
using the framework for quantum gravity adopted in namely that of Liouville 
strings in which classical conformal string backgrounds are allowed to exhibit 
quantum space-time fluctuations. These cause departures from criticality and con- 
formal symmetry that may be counterbalanced by introducing on the string world 
sheet a Liouville field with non-trivial dynamics. This Liouville field may be consid- 
ered as a dynamical renormalization scale, and conformal invariance may be restored 
by dressing operators with appropriate Liouville factors ||^. Since the quantum- 
gravitational space-time fluctuations drive the string super-critical, the Liouville 
field has negative metric, and we identify it with target time 0. 

Interactions of the Liouville field with the conventional string degrees of freedom 
are controlled by the Zamolodchikov flow of the world-sheet renormalization- 
group equations. These incorporate the conventional Hamiltonian structure and 
S-matrix description of the scattering of low-energy particles, but also encode infor- 



mation on the energy- dependent interactions of the probe with non-local quantum- 
gravitational degrees of freedom. These fail to decouple |0 from the low-energy 
matter in the presence of foamy singular metric fluctuations. Since a low-energy ob- 
server cannot detect such global modes by direct scattering experiments, integrating 
them out of the low-energy effective theory makes the low-energy system resemble 
an 'open' quantum system interacting with an unobserved 'environment' of these 
quantum-gravitational solitonic states. 

The outline of this paper is as follows. In Section 2, we review relevant features 
of this Liouville string framework for space-time foam, discussing in particular 
the formal basis for the appearance of stochastic fluctuations in theory space. Next, 
in section 3, we construct a specific realization of this non-critical string approach 
by considering quantum fluctuations of D-brane excitations in the vacuum. Then, 
in Section 4, we discuss the propagation of photons through this quantum D-brane 
foam. Section 5 shows how stochastic fluctuations occur in the context of D branes, 
and exhibits the expected diffusion effect Q Finally, in Section 6, we discuss how this 
phenomenon may be constrained or detected by observations of distant astrophysical 
objects. 



2 Relevant Aspects of Liouville String 

In this Section we review briefly features of Liouville that play roles in our subse- 
quent discussion. In this approach, one treats quantum fluctuations in the space- 
time background as deviations from a classical string background that is described 
by supercritical deformations of a conformal field theory on the world sheet. One 
may restore criticality by Liouville dressing of world-sheet model fields p. Specif- 
ically, consider a conformal a model, described by an action S* on the world sheet 
S, subject to non-conformal deformations g^Vid^a^ with Vi appropriate vertex 
operators: 

Sg = S* + J^g'V.d^a (1) 
The non-conformal nature of the couplings implies that their world-sheet renormal- 

ization-group (3 functions do not vanish. The following is the generic structure of 
such P functions, close to such a fixed point {g"^* = 0}: 

f3^ = (h,-2)g' + d^,g^g' + o{g'). (2) 

One cancels the deviation from criticality by world-sheet 'gravitational' dressing 
that corresponds to defining renormalized couplings in a curved space. To O^g"^), 
one has § |, |l3: 

X\t) = A'e"»* + ^ dfcA^A^te"'* + ©(A^), Q^ = c- 25 (3) 
^ ^ Q±2ai V ^ V / 

^We do not find any evidence to support the suggestion of birefringence |3|. 



where t is the zero mode of the Liouville field, is the central charge deficit which 
is > for the supercritical string case of interest here, and the ctj are gravitational 
anomalous dimensions: 



a,{a, + Q) =hi-2 for c > 25 (4) 



The supercriticality implies a Minkowskian signature for the Liouville field fTT 



enabling us to identify its zero mode with target time Q- After the renormal- 
ization (|^), the critical-string conformal invariance conditions corresponding to the 



vanishing of flat-space /5 functions are replaced by H, P, 10, 



A' + QX^ = -(3' fore > 25. (5) 

The minus sign in front of the fiat- world-sheet (3 functions reflects the supercriticality 
of the string. 



The propagation of non-relativistic light-particle modes was examined in where 
a modification of the quantum Liouville equation for the density matrix was found. 



as proposed in [0 . The propagation of massless probes was examined in [|[] , where 
it was found that the conventional relativistic energy-momentum dispersion relation 
is modified in non-critical Liouville strings, as a result of the interaction with the 
quantum-gravitational environment, as we now review briefly for the benefit of non- 
expert readers. 



The first step is to observe that, in the case of interest, the non-criticality of the a 
model describing the effective theory is induced by the operator-product-expansion 
coefficients c*- that express the interaction of a low-energy probe (latin indices) 

with quantum-gravitational modes (hatted latin indices). Hence, = to lowest 
order, and the Liouville anomalous- dimension coefficients are given simply by the 
magnitude of the spatial momentum of the massless probe: 

Oii = \k\ (6) 

We use this to rewrite (H) approximately, to order g"^, as |1|: 

where the shift Aa^ is given by: 

Aai ~ —c\g^ (8) 

^ We note for later use that in the gravitationaUy-dressed (LiouviUe) world-sheet theory, 
only the leading-order coefhcients in the /3 functions are renormalization-scheme independent. This 
implies that, when one identifies the Liouville field with target time, one loses general covariance 
in the foamy ground state. This should be thought of as a spontaneous breaking of the symmetry 
by quantum-gravitational fluctuations. 



(7) 



We next make the generic hypothesis, which is supported by some exphcit exam- 
ples [0, especially in the context of D branes [|14|, as discussed in more detail 



below, that 

-d. ~ iE/MQG, (9) 

where E is the energy scale of the low-energy probe, ^ = ±(9(1) and Mqg is a 
characteristic of gravitational interactions, possibly to be identified with the Planck 
scale Mp ~ 10^^ GeV. We infer from (H), (|^), (H) and (||) a modified dispersion 
relation 

i?= |fc| X {I + ^E/Mqc) (10) 
corresponding to an energy- (frequency-)dependent refractive index 77 = 1+0{E IMqg). 

It was observed in that such a quantum-gravitational effect could be probed 
by observations of distant astrophysical sources, such as GRBs, Active Galactic 
Nuclei (AGNs) or pulsars, and that these could place important constraints on 
models in which such a refractive index is suppressed only by a single power of 
1/Mqg, potentially with sensitivity to Mqg ~ Mp. Here we go further, remarking 



that quantum effects due to world-sheet topology fluctuations ||T^, which arise in 
the context of first-quantized strings from the summation over genera, may cause 
observable stochastic effects that lead to a diffusive spread in apparent velocities 
even for photons of fixed energy (frequency). This would imply the spreading of an 
initial pulse, and place limitations on the resolutions of experimental measurements 
of distant astrophysical sources. 

The basic motivation for this suggestion comes from the quantization of the target- 
space fields {g^} in string theory, which arises from the summation over world- 
sheet topologies. In our case, this procedure leads to quantum fluctuations in the 
cr-model couplings felt by the propagating (low-energy) string particle modes. 
The probability density V[g] in the space of cr-model theories is given in leading 
approximation by: 

V[g] ~ exp i-g'^g^ (11) 



where Gij ~< ViVj > is the Zamolodchikov metric in theory space for vertex 
operators Vi corresponding to the couplings g^. Formally, the non-trivial probability 



density (|Il|) arises |1T5|, as a result of the requiring the cancellation of modular 



infinities against renormalization-group divergences associated with singular con- 



figurations in the sum over genera ||15[. The width T in (|lT|) depends on some 
positive power of the string coupling constant gg, and the precise form is model- 
dependent [|16|, 0. The quantum uncertainties 6g^ in the couplings are found by 



diagonalizing the basis in theory space g'^, for which knowledge of the Zamolodchikov 
metric is essential. 



In the presence of fluctuations in the world-sheet topology of the string, which 
were not discussed in [0, ||, the shifts (§) associated with the refractive index of the 
vacuum will fluctuate as: 

|5(Aa,)|~|^c^5/| (12) 

leading in turn to stochastic quantum fluctuations in the refractive index in vacuo. 
This in turn leads to a characteristic diffusive spread in the arrival times of photons 
identical energies, to which we return later. Before that, we first compute the 
Zamolodchikov metric and the related width F (p!T| ) in a specific model for quantum- 
gravitational foam based on D particles. 



3 Fluctuations in i^-particle Foam 



We compute the quantum flucuations in ([IT]) in the particular case where the grav- 
itational degrees of freedom g'' are the collective coordinates and/or momenta of a 
system of N D particles |T^, |T^. As discussed in the hterature |TB|, |Ty|, |T^, [T^ such 



collective coordinates can be described by operators characterizing the recoil induced 
by the scattering of string matter off the D-particle background. The corresponding 
cr-model deformation is: 



Y^ix'') = isiY^'ise + f/fx")eee(x") (13) 

where the spatial coordinates of the D particle are identified with the couplings 
Yi and the Ut correspond to their Galilean recoil velocities. These coordinates are 
to be regarded as a-model couplings in the sense of the previous section. The 
parameter e — 0+ regulates the ambiguous value of 9(s) at s = 0, which ensures 
that the D-particle system starts moving only at the time x*^ = It is related to 
the world-sheet ultraviolet cutoff scale A (measured in units of the world-sheet size 
S) by = 2£^logA, where is is the fundamental string length. For finite e, the 
operators Yi, Ui each have an anomalous dimension — ||ep < ||T8|, and thus lead to 
a relevant deformation of S*. The corresponding renormalization-group equations 
T§: 

dY^^'^/dt = Uf, dUf/dt = 0, (14) 



are 



which are the Galilean equations of motion for the D particles, if we identify the 
time with the world-sheet scale: t = ig log A. 



The natural geometry on the moduli space Ai of deformed conformal field theories 
described by the above recoil operators (|13]) is given by the following Zamolodchikov 

■^Throughout this paper, we consider the ct- model 'renormaHzed' velocity ||l9|, |l^, |l^, which is 
related to the 'bare' velocity Usi by Ui — Usi/s- It is Ui which is an exactly-marginal deformation 
in a world-sheet renormalization group sense, and therefore corresponds to uniform i?-particle 
motion for times t > 0. 



metric 0^^= {V^^V^,) O: 



ab.cd 



+ ®U' + (U'U' + U'U') ® In}] + O (gf) 



db;ca 



(15) 



where In is the N x N identity matrix and we have renormahzed Qs to the time- 
independent couphng gs/\e\is- The canonical momentum P^^ of the D-particle sys- 
tem is given in the Schrodinger picture by the expectation value of —iS/dY""^ eval- 
uated in a cr model deformed by the operator V^f,, i.e., P*^ = {V^^) |T^: 



ab 



6 ^ 



oig. 



- ba 



(16) 



which coincides with the contravariant velocity P^^ 

In agreement with general arguments 0, we note that the moduli space dynamics 



can be derived from a Lagrangian of the form 



C 



'J_Yaby-<ij Y 
c) i ^ ab\cd j 



cd 



c 



(17) 



which coincides to leading order with the non-Abelian Born-Infeld effective 
action for the target-space D-particle dynamics: 

-Cnbi = T— Tr Sym /det [7]mn In + ^sdl Fmn] (18) 



s9s 



where the trace Tr is taken over U (N) group indices, 

Sym(Mi,...,M„) = ^ E ^-i " " " 



(19) 



is the symmetrized matrix product, and the components of the dimensionally- 
reduced field-strength tensor are given by Pqi = js^i and Fij = ^[Yi,Yj]. In the 
Abelian reduction to the case of a single D particle, the Lagrangian ( |18|) reduces to 
the usual one describing the free relativistic motion of a massive particle. The lead- 
ing order P^ term in the expansion of (|18D is just the usual Yang-Mills Lagrangian. 



The formalism described in is a non-trivial application of the theory of Liou- 
ville string and logarithmic operators. We note, moreover, that the derivation of the 
Lagrangian (|l8|) was made 'off-shell', in other words we have compared generalized 
momenta in theory space with those derived from the Born-Infeld lagrangian. The 
equivalence between the two formalisms extends beyond the equations of motion, 
which are the conformal-invariance conditions of the a model. This is a central 
aspect of the recoil approach: there are deviations from the usual equations of 
motion for low-energy modes, because conformal invariance is violated by the re- 



coil process |]T8|, |T9[ , and the Zamolodchikov world-sheet renormalization-group flow 
provides an off-shell treatment of this recoil problem. 



The leading contributions to the quantum fluctuations in Ai arise from pinched 
annulus diagrams in the summation over world-sheet genera of the a model . 
Symbolically, these lead to contributions in the genus expansion of the form 



•d + io + * + 



(20) 



consisting of thin tubes of width 5 — (world-sheet wormholes) attached to the 
world-sheet surface S. The attachment of each tube corresponds to inserting a 
bifocal pair V^f^{s)V^^{s') on the boundary dT,, with interaction strength g"^, and 
computing the string propagator along the thin tubes. There are modular diver- 
gences of the form log 6, which should be identified with world-sheet divergences at 
lower genera 0, so we set 

\ogS = 2g^JogA=l^g]e-' (21) 



The exponent r/ 7^ in this Fischler-Susskind-like |^T| relation, which allows one to 
cancel logarithmic modular divergences by relating the strip widths S to the world- 
sheet ultraviolet scale A, arises because the relation (pT]) is induced by the string 



loop expansion. As we have argued in |jT^, the case 1] ^ allows direct compar- 



ison of our results with other results in the string literature, based on alternative 



approaches |22]. However, we do not determine t] here, since we only consider string 
interactions between D branes. Considering brane exchanges between the system of 
D branes, may enable the value of 77 to be fixed (see also below), but this point is 
not important for the present analysis. 

One effect of the dilute gas of world-sheet wormholes is to exponentiate the bilocal 



operator, leading to a change in the a- model action |T5|, |T6[. This contribution can 
be cast into the form of a local action by rewriting it as a Gaussian functional 
integral over wormhole parameters as described in the previous section, and we 
arrive finally at [17]: 



E 

genera 



Zj,[Y] ^(1 Dp ^f^'G^,^..pfm^^lai^o^5 ^^q^. y ^ A (22) 
\jm 1 



We see from (^2]) that the effect of this resummation over pinched genera is to in- 
duce quantum fluctuations of the solitonic background, giving a statistical Gaussian 
spread to the D-particle couplings. Note that the width of the Gaussian distribution 
in (p2D, which we identify as the wave function of the system of D particles 0, is 
time-independent, and represents not the spread in time of a wave packet on A^, 
but rather the true quantum fluctuations of the D-particle coordinates. 



The corresponding spatial uncertainties can be found by diagonalizing the Zamolod- 
chikov metric (|T5|), as was done in [|l^ and will not be repeated here. For the case 



of a single D particle: a = b, one arrives at the uncertainties: 



AX! 



= AY' 



sils 



1 + 



87r3 



2 



7?/2 



(23) 



for the individual D-particle coordinates, 
coincides with the standard string smearing 



For rj 

m 



= 0, the minimal length in (|23|) 
whereas for r/ = | it matches the 
form of ip which arises from the kinematical properties of D particles ||22|. A value 
7] 7^ is more natural, because the modular divergences should be small for weakly- 
interacting strings. We note that the uncertainty ( p3D is time-independent 



which is important for experimental tests of the phenomenon, as we discuss below. 
The coordinate uncertainties for a ^ b are responsible for the emergence of a true 
non- commutative structure of quantum space-time, and represent the genuine non- 
Abelian characteristics of multi-D-particle dynamics, but we are not concerned with 
such a case here. 



4 Propagation of Photons in a D-Particle Foam 
Background 

After this general discussion, we now discuss the propagation of photons in the 
background of multiple D particles, taking into account the recoil of the latter due 
to the scattering of the photons. This is immediate in the context of [|l^: one 
simply adds to the argument of the determinant (plSf ) an Abelian U{1) field strength 
fhiN = dMdN — Qncim, where a*/ denotes the U{1) electromagnetic potential of 
the photon fields. The corresponding Born-Infeld effective action is considered in 
four-dimensional space time. We adopt the conventional view point that photons 
interacting with D particles may be represented by adding to the cr-model action So- 
on the (open) world sheet an electromagnetic potential background as a boundary 
term [||]: 

S^3- I dTaM{X)drX^' (24) 
c jas 

in a standard (Neumann) open-string notation ^ 



The non-relativistic heavy D-particle []r^ corresponds to a time-dependent back- 
ground ([13|) , whilst the photon propagator depends on the full four-dimensional 
space time. It is important to note, though, that in our approach the coupled 
Z^-particle-photon system is out of equilibrium, due to the recoil process and its as- 
sociated distortion of the surrounding space time. This is reflected in the fact that 



^The target-space Born-Infeld action, including photons interacting with D particles, may al- 
ternatively be considered as the target-space action of a three brane, i.e. a solitonic object in 
string theory with four coordinates obeying Neumann boundary conditions and the remainder 
Dirichlet boundary conditions [|l8j . It is possible that the above Neumann picture can be obtained 
from the Dirichlet one by a world-sheet T-duality transformation, but there are problems with this 
transformation at a quantum level [E5[ . We restrict ourselves here to the Neumann string picture. 



the resulting backgrounds do not satisfy the classical equations of motion, as was 
mentioned earlier in the general context of the departure from conformal symme- 
try. Moreover |TB[, the recoil curves the surrounding space time, since it induces - 
via Liouville dressing and the identification p, |13| of the Liouville mode with the 
target time t - graviton excitations for the string, which in a cr-model framework 
correspond to graviton backgrounds with non-trivial off-diagonal elements: 



(25) 



where f/j is the velocity of the recoiling D particle. 



The element (p5|) is obviously not Lorentz covariant, reflecting the spontaneous 
breaking of this symmetry by the ground state of the string. The splitting between 
the quantum-gravitational 'medium' and the propagating particle subsystem is not 
possible in a Lorentz-invariant way, and there is no formal reason that this should 
be so. On the contrary, spontaneous violation of Lorentz symmetry is generic in 
Liouville strings |Tl|, ^. Thus we may expect the interaction of such graviton modes 
with the photons to lead to a modification of the photon dispersion relation, anal- 
ogous to non-Lorentz-covariant, e.g., thermal, effects in conventional media. A key 
difference between quantum-gravitational effects and those in conventional media 
is that the former increase with the energy of a probe, whilst the latter attenuate 
with increasing energy. 



Gravitational interactions may be incorporated in the Born - Infeld lagrangian 
|8D by replacing t/mn ^ Gmn and using the non-covariant expression (pSf ) when 



calculating amplitudes. Thus the complete off-shell Born - Infeld Lagrangian for the 
interaction of photonic matter with D particles is: 



C 



NBI 



— Tr Sym 



s9s 



\ 



det 

M,N 



Gmn In + — fivmio) In + ^IqI Fmn 
c 



(26) 



which is the basis for our subsequent discussion. 



We first re-examine the possible order of magnitude of the frequency-dependent 
refarctive index (H) induced by such effects, concentrating first on the case of Abelian 
(single D-particle) defects. The terms in the effective action ([T8| ) that are relevant 
for the modification (^ of the photon dispersion relation arise from the three-point 
function terms in (^ that involve two photon excitations and one induced-graviton 
excitation (^). The appropriate term in a derivative expansion of the Born- Infeld 
action is: 

j^photons, graviton ^ f j^[a)G N AM ~ fij{a)UjfoieH (27) 

where latin indices are spatial, and we used (pSf). Terms that are similar in order 
of magnitude, but not in tensorial structure, can also be obtained by combining the 



fliN terms with the determinant \J —det{GMN) of the target-space metric f\ An 
appropriate order-of-magnitude estimate of the recoil velocity is: 



Uir-. g^\k\/Mp (28) 

reflecting energy- momentum conservation in the scattering process IB, 0, and 
assuming that the momentum of the recoiling D particle is of the same order as 
the energy-momentum of the photon. Working at times t after the recoil that are 



sufficiently large, we may assume |jT8|, |T9|, |T6|, 0] that 

"2 - t (29) 



on average. Using now for the graviton a linear approximation about flat Minkowski 
space Gmn ~ VMN + huN, which allows one to consider a simple Fourier momentum- 
space decomposition of the pertinent amplitudes, we conclude that the effective 
photon-graviton interactions in (^) are of order 0{gsE\k\'^/Mp). Such estimates 
apply to the string amplitudes c\^g'' in the shift (|^) in the single-defect case, which 
therefore becomes, to leading order in a low-energy approximation: 

Aa. = 0{-gs^)) (30) 

Using (|^), we see immediately that this leads to a modifled dispersion relation for 
the photon: 

(31) 

The refractive index is then determined from the photon group velocity 0: 

The sign of the refractive index is determined by the fact that the Born-Infeld action 
underlying the above analysis prevents superluminal propagation. This sign was not 
determined in our previous discussion @]. As discussed there, such a variation in 
the velocity of light will cause a spread in the arrival times of pulses of photons, 
according to their energy p[: 

CLE , , 

A* ~ ^ (33) 

The possibility of observing experimentally such a shift using distant astrophysical 
sources appears conceivable, as dicussed in and reviewed in the last section of 
this paper. 



'^It is clear that Lorentz-invariant interactions of powers of p terms with generic Tr(F^") terms 
in the derivative expansion of ( p6|) do not affect the photon dispersion relation. It is only non- 
covariant terms, such as the the induced-graviton-photon interactions considered above, that do 
so. 



5 Stochastic Fluctuations in the Apparent Veloc- 
ity of Light 

We now extend the above discussion to include an estimate of the stochastic fluc- 
tuations in the apparent velocity of light, and hence the refractive index, due to 
the summation over higher-genus world-sheet topologies. The general theory of this 
summation in the context of Liouville string, discussed in section 2, leads to stochas- 
tic fluctuations in the collective coordinates of the D particles discussed in section 



3 0: 

SY' > e^g]/', (34) 

The conventional Heisenberg uncertainty relation between the coordinates and 
the corresponding canonical momenta Pj has been shown to take the following form 
for D particles, after the summation over genera ||17|| : 

6Y'6Pj > 2g]+'^l^5] (35) 



Saturating the bound (|34D in (|35D , we obtain the following estimate of the uncer- 



tainty in the associated collective canonical momenta of the D particles: 

m ^ (36) 

We note that this estimate of the uncertainty in the collective momentum is inde- 
pendent of the exponent ?7, which, as already noted, would be required to take the 
value = 2/3 in order to match results on D particles in conformal string theory. 



In our interpretation, the fluctuations arising from the summation over genera 
lead to a statistical superposition of theories with different values of the couplings 
gK This 'stochastic environment' is characterized by the Gaussian form ( |TT1) of the 
associated probability distribution in the space of cr-model backgrounds. A pho- 
ton of given energy E propagating through the fluctuating quantum-gravitational 
medium is subject to stochastic fluctuations in its velocity, which may be obtained 
from Liouville interactions of the form ([l^) by endowing the generalized cr-model 
coordinates {g^} with the fluctuations (^) found for the collective coordinates Y^ 
of the D-particle foam. 



To see how this effect indeed induces stochastic fluctuations in the refractive index 
of the photon, one first calculates the amplitudes c^-j^g^ appearing in ([T^) using the 
Born-Infeld Lagrangian (|26|) , calculated on world sheets with disc topology. Using 
(|29D, one sees that the leading terms in a derivative expansion are the ones given in 
(|27D. The summation over world-sheet genera leads to stochastic fluctuations in Ui 
(|28|), which are given by (^): 



(37) 



These stochastic fluctuations in Ui generate fluctuations in the overall proportion- 
ality coefficient of the Maxwell terms in the effective action for the photon: 

^ ^photons, graviton ^ J^^^(q) {5Gna) f AM ~ fij{o) (SUj) foieH (38) 

whose Fourier transform is 0{2g^E\k\/Mp). This does not itself affect the velocity 
of the photon, but simply renormalizes the energy scale. 



However, when one proceeds to higher orders, one picks up contributions that 
lead to stochastic fluctuations in the refractive index. To see this, we concentrate 
on terms of quadratic order in Ui, stemming from terms in a derivative expansion of 
the Born-Infeld lagrangian ( PUj) that are of the generic form /I^^TrF^^, where the 
non-Abehan field strength F^f^ is calculated by considering Y"''' (|l^) formally as a 
'gauge potential'. For long times of order (pQ]), the only non- vanishing components 
of Fab are Foi ~ Ui which lead to stochastic fluctuations in the quadratic terms of 
order 

fllN {SU^) = ^9lfMNU^ (39) 

Thus, in order to discuss the leading-order effects induced by the stochastic fluc- 
tuations of our D-particle foam, as implied by the summation over world-sheet 
topologies, we should consider corrections to the Liouville-dressed couplings that go 
beyond quadratic order in the a-model couplings {g^}- Such corrections have been 
studied in ||10|], and are not given here explicitly. It is sufficient for our purposes 
to point out that some are proportional to four-point amplitudes, c*;,^ divided by 
'energy denominators' aj+at'- c^ji^J [ai+Uj), in the limit of vanishingly small central- 
charge deficits Q that are of interest to us, whereas others are of the form c*„jC^ 
divided by terms of order af. The three- and four-point amplitudes are computed 
from the Born-Infeld lagrangian (|26|) above 0. 

Thus, for long times of order (PP|), and using ( P^ ) for Ui, one finds the following 
estimate for the corresponding stochastic fluctuations in (|]): 

MAa.) = ^ (40) 

Correspondingly, we find fluctuations in the velocity of light in the quantum-gravitational 
medium of order 5c ~ Sg^E/Mp, motivating the following parametrization in the 
stochastic spread in photon arrival times: 

{6 At) ~ ^ (41) 



where Aqg ~ Mp/Sg^. We emphasize that, in contrast to the variation (33) in 
the refractive index, which refers to photons of different energy, the fluctuation (|4l| ) 
characterizes the statistical spread in the velocity of photons of the same energy. 

^We again reminder the reader that, since such higher-order corrections are world-sheet 
renormahzation-scheme dependent ]To| , this reflects the (spontaneous) breaking of general co- 
variance by our foamy ground-state when the Liouville field is identified with target time M. 



6 Observational tests 



The most important signatures of the refractive index and the stochastic fluctuation 
in the velocity of light that we find are that they increase linearly with the photon 
energy (frequency). This means that they can, in principle, easily be distinguished 
from more conventional medium effects, that attenuate with increasing energy. The 
effects scale inversely with some quantum-gravitational scale characteristic of strings 
and D branes. We are not in a position to estimate it numerically, but we expect 
it to be within a few orders of magnitude of Mp ~ 10^^ GeV. In standard string 
theories one has 77 = 2/3 and gg/4:iT ~ 1/20, but the latter may well be modified 
in a more realistic theory. In principle, one could even envisage using upper limits 
on (measurements of) the rate of broadening of a radiation spike of definite energy 
(frequency) to constrain (measure) gs- 



We conclude this paper by mentioning some possible observational tests of these 
ideas. As has been emphasized previously 0, the figure of merit for constraining 
the possibility of an energy- (frequency-)dependent refractive index in vacuo is the 
combination L x AE/At, where L is the distance of a source, AE is the range 
of photon energies studied, and At is the observational sensitivity to differences in 
light-travel times. The latter is limited by the durations of pulses produced by the 
source as well as by the resolution of the detector. The corresponding figure of merit 
for testing the new possibility advanced in this paper, namely a stochastic spread in 
light-travel times for different photons with the same energy E, is simply L x E / At. 
In practice, when comparing photons of different energies, AE is often dominated 
by the highest photon energy E that is measured, so that AE ~ i?, and the two 
figures of merit are essentially equivalent. 



Astrophysical sources offer the largest figures of merit, and the most promising 
that have been considered include GRBs 0, AGNs ||2^ and pulsars [Q. These 
have all been used already to constrain the refractive- index parameter Mqg, and 
offer similar sensitivities to the stochastic-spread parameter Aqg- As mentioned 
earlier in this paper, the D-brane Born-Infeld analysis indicates that higher-energy 
photons should be retarded relative to lower-energy photons, rather than advanced, 
and their arrival times should be more spread out. In the Table below, we list some 
of the sources that have been considered, and the sensitivities (limits) that have 
been obtained. For completeness, we have also indicated the sensitivity that might 
be obtainable from a detailed analysis of the recent GRB 990123. We see from the 
Table that Mqg cannot be much smaller than the Planck scale Mp ~ 10^^ GeV, 
and that some of these astrophysical sources may already providing sensitivities to 
Mqg, Aqg ~ 10^^ GeV. This provides additional fundamental-physics motivation 
for such 7-ray observatories as AMS and GLAST PD|. 



Other probes of the signatures of quantum gravity that might be provided by 
the unorthodox photon propagation proposed here might be possible using labo- 



ratory experiments, for example those testing quantum optics and searching for 
gravitational waves, but we do not explore these possibilities further in this paper [|. 
However, we think that the discussion given here demonstrates amply the possibility 
that at least some quantum-gravity ideas may be accessible to experimental test, 
and need not remain for ever in the realm of mathematical speculation. 



Table: Observational Sensitivities and Limits on Mqg, Aqg 



Source 


Distance 


E 


At 


Sensitivity (Limit) 


GRB 920229 


3000 Mpc (?) 


200 keV 


10-2 s 


10^6 GeV (?) 


GRB 980425 


40 Mpc 


1.8 MeV 


10-3 S (?) 


10^6 GeV (?) 


GRB 920925c § 


40 Mpc (?) 


200 TeV (?) 


200 s 


10^9 GeV (?) 


Mrk421 |7| 


100 Mpc 


2 TeV 


280 s 


> 4 X 10^6 GeV 


Crab pulsar p8|| 


2.2 kpc 


2 GeV 


0.35 ms 


> 1.8 X 10^^ GeV 


GRB 990123 


5000 Mpc 


4 MeV 


ls(?) 


3 X 10^^ GeV (?) 



The question marks in the Table indicate uncertain inputs. Hard limits are indicated 
by inequality signs. 
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